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ABSTRACT 

In a preceding paper, using Eggleton's evolution code we have discussed the structure 
and evolution of low-mass W UMa type contact binaries without angular momentum 
loss (AML). The models exhibit cyclic behavior about a state of marginal contact on 
a thermal time-scale. Part of the time of each cycle is spent in contact and part in 
a semi-detached state. According to observations, W UMa systems suffer AML. We 
present the models of low-mass contact binaries with AML due to gravitational wave 
radiation (GR) or magnetic stellar wind (MSW) are presented. We find that gravita- 
tional radiation cannot prevent the cyclic evolution of W UMa systems, and the effect 
of gravitational radiation on the cyclic behavior of contact binary evolution is almost 
negligible. We also find that the most likely AML mechanism for W UMa systems is 
magnetic braking, and that magnetic braking effects can increase the period of the 
cyclic evolution, and shorten the fraction of the time spent in the poor thermal contact 
state exhibiting EB light curve. If W UMa stars do not undergo cyclic evolution, and 
their angular momentum loss is caused simultaneously by MSW of both components, 
we find that the value of the parameter. A, should be taken a larger value in compari- 
son with those derived from observations of single stars. This indicates that the AML 
efficiency in W UMa systems may be lowered in comparison with non-contact stars 
because of the less mass contained in the convective envelopes of the components in 
W UMa systems. If W UMa systems lose their angular momentum at a constant rate. 
An angular momentum rate of « 1.6 x 10"^yr~^ can prevent the cyclic behaviour 
of the model, and the model can keep in good contact with an essentially constant 
depth of contact. 
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1 INTRODUCTION 

In a preceding paper (Li et al. 2004, hereafter Paper I), the 
structure and evolution of low-mass W UMa systems have 
been dis cussed by us using the stella r evolution code devel- 
oped bvlEggletonI Jl97ll.ll97lll973l) . with updated physics 
JPols et alJll995^ . Our models for low-mass W UMa sys- 
tems were based on the conservative assumption, namely 
that the total mass of the system and its orbital angular 
momentum remain preserved. As a result, our models ex- 
hibit cyclic behavior about a state of marginal contact with 
a period of about 10^ yr. Part of the cycle is spent in con- 
tact and part in the semi-detached state. In any cycle, a 
larger temperature difference (ATch > 300K) between the 
two components occurs in a part of a cycle (about 30-35 per 
cent time of a cycle). So most previous investigators had 
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thought that this requires there to be as many short-period 
binaries with EB light curves as with EW light curves, and 
that the TRO models for W UMa type systems cannot ex- 
plain the light curve parad o x according to the investiga- 
tion of lLucv fc WilsonI ill979r i. lRucinskil ||2002^ finds 13 EWs 
and 5 EEs [and 14 ellipsoidal (ELL) variables, which have 
too small an amplitude to be classified as EWs or EEs)]. 
It is reasonable to identify the EWs as contact binaries, 
and the EEs as semi-detached. The ratio of 5/13 is not 
much out of line with TRO theory l! Lucylll976l: iFlannervl 
1976; Ro bertson fc Eggletonl ll977'l indicating that zero-age 
contact systems are not thermally stable, but suffer cyclic, 
thermally unstable mass transfer. 

The TRO models for W UMa systems (Rahunen 1981, 
see also Paper I) has shown that the predicted ratio of 
the timescales characterizing the phases in which the os- 
cillating systems will exhibit W UMa-type (EW) and /3 
Lyrae-type (EB) light curves is tew/teb < 2 (in the mass 
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ratio range 0.5 < g < 0.7). But the observations give 
tew/teb ~ 5 to 6, suggesting that W UMa systems proba- 
bly suffer angular momentum loss. It is well known that loss 
of angular momentum from the binaries consisting of two 
point masses via the emission of gravitational wave radia- 
tion (GR) causes the separation between the components 
to decrease. The possible importance of gravitational ra- 
diation for short period bin aries was first pointed out by 
iKraft. Mathews, fc Greenstein 0,962,) • In the evo l ution of W 
UMa systems it has been d iscuss ed bv lFaulkneil Jl97J) and 
IPaczviisd l(l963). iHuand il966l) suggested that magnetic 
torques could bring together the separate components of a 
detached binary. This is Schatzman's (1962) and Mestel's 
(1968) mechanism of magnetic braking. If the rapid rotation 
leads to strong dynamos in the convective zones of the solar 
stars in W UMa systems, effective magnetic braking may 
result and the angular momentum loss by magnetic stel- 
lar wind (MSW) may control the contact binary evolution 
l|^n]r^3[l97i). Meanwhile, there is increasing evidence 
of strong magnetic activity in short-period binaries, includ- 
ing W UMa stars (spots, flares, strong chromospheres and 
coronae, etc.). This suggests indirectly tha t W UMa systems 
suffer angular momentum loss via MSW. lOkamoto fc Satol 
l|l970fl have considered the forma tion of W UMa stars by 
magnetic braking in more detail. IVan't Veed lll979l 1198(1) 
has further developed this idea. The angular moment um los s 
in contac t binar ie s has also been treated bv ,Mosj 1197^ , 
IWebbinkI (Il976l) iRobertson fc EggletonI Jl977ft . iRahuned 
||1981^ and TvilhuHl982^ ! 

The observational data and theory of cont act binaries 
have been rev i ewed in thr ee extensive p apers bv lMochnackil 
Jl98ll) . IVilhl Jl98ll) and lSmithI il984t) . Currently the pre- 
vailing opinion is that the AML phenomenon is either the 
main property of or an ingredient in the proper explanation 
of the existenc e and evolution of contact binar ies of late 
spectral types dVan't Veejll979l : iRuciriskil [1983) . There is 
ample evidence that the W UMa-type systems posses all the 
chromospheric and coronal phenomena normally (i.e. in the 
solar analogue) related to the existence of magnetic fields. 
Robertson & EgEloton ( 1977^) have found that a steady loss 
of angular momentum drive s the system towards smaller 
mass ratios. iRahunenI il98 jl has found that it is possible 
to keep zero-age models in good thermal contact if the or- 
bital angular momentum of the binary is allowed to decrease 
and the timescale of the required angular momentum loss is 
about 5 X 10* yr. 

In this paper, we restrict our attention to the general na- 
ture of the effects of angular momentum loss (i.e. J < 0) and 
to the influence of the matter accreted by the secondary on 
its evolution in a semi-detached state and present models of 
low-mass W UMa systems with angular momentum loss via 
gravitational radiation or magnetic braking. We find that 
the infiuence of gravitational radiation on the cyclic evo- 
lution of contact binaries is almost negligible because the 
gravitational radiation life-time is long compared with the 
thermal timescale of the primary of our models. Systemic 
angular momentum loss via gravitational radiation only ac- 
celerates the evolution towards more extreme mass ratios, 
rather than leading to orbital collapse. The most likely an- 
gular momentum loss mechanism for W UMa systems is 
magnetic braking. It can increase the period of the cyclic 
evolution and shorten the fraction of the time spent in the 



state exhibiting EB light curves, and a suitable angular mo- 
mentum loss rate can prevent the cyclic behaviour of the 
model. Meanwhile, if we consider an energy source at sec- 
ondary's atmosphere provided by the accreting matter from 
the primary in semi-detached evolution, and find that the 
expansion in radius of the secondary is hastened by the 
accreting matter which should have a higher entropy than 
the original star, it can shorten the time spent in the state 
exhibiting EB light curve. We compare the convective en- 
velopes of a semi-detached model with those of a contact 
one which has the same total mass and mass ratio as the 
semi-detached one. We find that the convective envelope of 
the secondary in the contact system is much thinner than 
that in the semi-detached one and the convective envelope 
of the primary in the contact system is slightly thicker than 
that in semi-detached one. But the mass contained in the 
convective envelopes of two components of contact system, 
which is related to the dynamo action, is less than that con- 
tained in the envelopes of the semi-detached one. This is 
caused by the energy transfer between the two components 
in contact systems. As a result, W UMa systems show the 
lower activity in comparison with the non-contact systems. 



2 CONTACT CONDITION AND LUMINOSITY 
TRANSFER 

2.1 Contact condition 

In Paper I, based on Roche geometry, we have given a radius 
grid of contact binaries. The radius grid gives the relative 
radii of both components of contact binaries with different 
mass ratios (0.02, 0.04,. ..,1.0) and different depth of contact 
(0.0, 0.025, 0.05,.. .,1.0) and it has been used as the surface 
condition of contact binaries to ensure the surfaces of both 
components lie on the same equipotential by interpolation. 
Here we use the same contact condition as Paper I. If the 
mass is transferred from star 2 to star 1, the rate of mass 
transfer is expressed as. 



— = -CMax[0,(ln — , 
in which 



^20 = 



(1) 



(2) 



The explanations of Eqs. (1) and (2) are as in Paper I. If 
mass is transferred from star 1 to star 2, the rate of mass 
transfer can be derived by the same way. 

2.2 Luminosity transfer 

2.2.1 Luminosity transfer due to mass transfer 

In the semi-detached phase, the mass above the critical 
lobe, together with the energy (including gravitational en- 
ergy, heat energy, and radiative energy) in the transferred 
mass, is transferred to the secondary. The effect of heating 
the photosphere of a low-mass ZAMS stars (with a deep 
convective envelope or fully convective), by the kinetic en- 
ergy of the infalling matter, is very important for their re- 
sponse JPrialnik fc Livid Il985^ . Since the radiative energy 
is so small that it can be neglected relative to gravitational 
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energy and heat energy, the accretion luminosity due to ac- 
creting mass is 



varies through the cycle and goes to zero with the depth of 
contact. Thus we take: 



iacc = I3{i>l - '02 + Si - S2)M2 



(3) 



where tpi and Si are, respectively, the gravitational poten- 
tial and specific entropy at the surface of the component 
i, M2 is the rate of mass accretion of the secondary and /3 
is a coefficient which describes the efficien cy of heating o f 
the photosphere by the shock wave region. iKahleil ll2002al) 
took the parameter /? to be 1.0, i.e. he thought that the 
bulk kinetic energy gained by falling through the potential 
is completely transformed into thermal energy and it, to- 
gether with the thermal energy carried by infalling matter, 
is regarded as an energy source in the outermost layer of 
the gainer without energy loss. However, some of this en- 
ergy may be dissipated dynamically in a shock, but the dis- 
sipated rati o of this energy is unkno wn, so we used a value 
of /3 = 0.5 iSarna fc Fedoroval[l989l) . This suggests that 50 
per cent of the transferred energy is lost from the system. 
Since the accretion luminosity is caused by energy sources in 
the gainer's outermost layers, we modify the surface bound- 
ary condition of the secondary to model the influence of the 
accretion luminosity on the evolution of the secondary, i.e. 
the surface temperature of the secondary may be approxi- 
mated by the formula, Lin + ^acc = inR^aT^g, where Lin 
is the luminosity coming to the photosphere from the stel- 
lar interior, R is the radius of the secondary, and a is the 
Stefan-Boltzmann constant. 



2.2.2 Luminosity transfer due to circulation currents in 
the common envelope 

IStruvd (^4^ first recognized the unusual mass-luminosity 
relationship of the secondary components of W Ursae Ma- 
joris systems, i.e. the two components have identical tem- 
peratures and hence luminosities scaling as L oc R^ (oc M 
for Roche geometry) in spite of differing masses. This Sug- 
gests that it might b e causally related to a possible com- 
mon envelope. lOsakil (HqG^) noted that von Zeipel's theo- 
rem would require the observed approximate constancy of 
radiative flux over the surface of a system with a radiative 
common envelope in order to maintain hydrostatic equilib- 
rium. The fact that most W UMa systems app ear to have 
convective envelopes, however, led ILucvI il968l) to propose 
that the anomalous mass-luminosity relationship possessed 
by W UMa binaries can be explained if energy transfer is 
assumed to take place between the two components of the 
systems well below the atmosphere. 

A luminosity transfer by circulation in the common en- 
velope is applied in the outermost layers which are close to 
the surface of each component. The amount of the luminos- 
ity transfer by circulation m ust satisfy the requirement of W 
UMa systems' observations (iRobertson fc EggletonllOT^ . it 
can be written as 



Li - ALo L2 + ALo 
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where L\^2 are the core luminosities (including nuclear and 
thermal luminosities) of both components, and mi, 2 the 
masses of the components. Since transfer is not fully effi- 
cient at all phases, an arbitrary factor / is introduced, which 



AL = /-ALo, (0</<l). 



(5) 



where / is the efficient factor of energy transfer. We take 
/ = Min[l,Q(d2-l)] (6) 



in which 
d = Max[l,Min 
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Rc-ritl ficrit2 

where ri,2 are the radii of both stars, -Rcriti,2 the Roche 
critical radii of both stars. The parameter a is expected 
to be moderately large, so that heat transfer becomes fully 
efficient for stellar radii exceeding the Roche radii by some 
standard small amount. In this work, we take a — 15. 

The luminosity transferred by circula tion cur r ents fr om 
the primary to the secondary adopted bv iKahleil (jioO^) is 

AL = / o-cx,2dm2 = - / acx,idmi, (8) 
Jo Jo 

where (Tcx,i is the source (when positive) or sink (when nega- 
tive) of energy per unit of mass caused by interaction of the 
components. Although the luminosity increment is applied 
in the adiabatic portion of the envelope of each star by most 
of the pre vious investigator s, according to the energy trans- 
fer model l|R,obertsonlll98n^ there is no essential distinction 
between convective and radiative envelopes. Meanwhile, Li 
et al. (Paper I) argued that the energy transfer may take 
place in the outermost layers for low-mass W-subtype sys- 
tems and the transfer takes place in the base of the common 
envelope for A-subtype systems. It follows that convection is 
by no means essential to heat transport in contact envelopes, 
although it may well have an important influence. There- 
fore, we assume that the luminosity transfer takes place in 
the outermost layers of the common envelope for low-mass 
W UMa systems, and the luminosity increment is applied 
in the 10 meshpoints which are close to the surfaces of both 
components as Paper I. 



3 GRAVITATIONAL RADIATION AND 
MAGNETIC BRAKING VIA MSW 

It is well known that loss of angular momentum from the 
binary systems consisting of two point masses due to the 
emission of gravitational radiation ca uses the separation 
between the components to decr ease ("Landau fc Lifschitzl 
■ 1962l) . .Faulkner (1971) and iTaa m. Flann crv, fc Fa ulknej 
1198(11 showed that gravitational radiation, by removing 
orbital angular momentum, could stimulate mass transfer 
from the main sequence red component at a rate which 
causes the binary to evolve more rapidly than nuclear cvolu- 
tion w ould dictate. Faulkner (1976) and Cliau & Lautcrborn| 
(Il977fl have both found very significant effects of J [via grav- 
itational radiation (GR) in their computa tions! on the evo- 
lution of low-m£iss close binary systems. ICharJ (Il978f) has 
considered the effects of loss of orbital angular momentum 
via gravitational radiation on the period changes and the 
implications for mass-transfer rate determinations. 

We adopt the standard Einstein theory which predicts 
that the rates of change of semi- major axis and orbital angu - 
lar momentum are, respectively, (iLandau fc Lifschitzlll962ll . 
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and 

din J _ 32G^ Ml Ma (Ml + Ma; 



(9) 



(10) 



At 5c5 yl4 

where J is the orbital angular momentum, A is the semima- 
jor axis of the system, Mi_a are the masses of both compo- 
nents, G is the gravitational constant, c is the light velocity, 
and t is the time. The rate of angular momentum lo ss associ- 
ated with gravitational radiation is also given by (iWebbinkI 
Il975l). 



dlnJ 



= -8.1 10 



_ 10 Ml Ma (Ml -I- Ma 



(11) 



At A^ 
where J is the orbital angular momentum in erg • s. Mi, a 
are the masses of both components in units of solar masses 
Mq, a is the semimajor of the system in unit of solar radii 
Rq, t is the time in years. 

Magnetic braking of the rotation of the cool stars in the 
binaries is an attractive mechanism. The idea that in stars 
with convective envelopes a dynamo is operating, which has 
efficiency that decreases with decreasing rotation rate, is 
consistent with observational data which indicate that the 
chromospheric (i.e. Ca II H and K) and coronal emissions 
are strongest in rapidly rotating stars l|Zwaan 198]|). The 
magnetic braking is presumably due to the stellar wind 
which streams out along the coronal ma gnetic field lines, 
which are rooted in the convection zone l|SchatzmarJll962L 
^6^. The loss of angular momentum per unit mass in 
the stellar wind is very large since the outflowing matter 
is forced by magnetic field to corotate with the star out 
to large distances. The magnetic stellar wind (MSW) as 
a possible driving force for the evo l ution of low-mass bi- 
naries was proposed first by iHuand l)l96ffl . Observational 
support for the action of such a wi nd in low-mass spectro - 
scopic binaries has been provided bv lKraicheva et all l)l97?jl . 
Theoretical studies of the effect of MSW on the evolution 
of close binaries consisting of a low-mass component with 
a convective envelope filling its Roche lobe and a com- 
pact accreting component that is a de generate dwarf or 
a neu tr on star have been conducted bvlVerbunt &: ZwanrJ 
IT^. lR,a,pna,nort. Verhunt .losj il9S.^). iTaaml jl 98.-j^. 



lR,a,nna,nort. Verhunt .los^ il9S.^) iTaamUl 98.-f . 
ml (11984). ISpruit fc Ritteii (Il98,^ ). iTutukovl 



IPattersor 

and lSarna &: Fedoroval |1989S), among others. 

The most likely angular momentum loss mechanism for 
W UMa systems is magnetic braking. The observational evi- 
dence and theoretical reasons have been presented for believ- 
ing that low-mass W UMa binaries can lose their orbital an- 
gular momentum efficiently as a consequence of a MSW from 
one or both components. The Emstem soft X-ray observa- 
tions and lUE ult r aviolet observations ([yaianalliflSft IVilhxJ 
1198,?!: lEatonlll98,'j: iRuciriski fc Vilhulll98,'j) showed that W 
UMa systems are strong sources. This suggests surface ac- 
tivity of the kind observed on our Sun, and so the pres- 
ence of magnetic field. A magnetic stellar wind would cause 
braking. In the absence of an explicit the ory for the angular 
momentum loss by MSW, we employ a Skum anichI (|l972f l 
law, Kot ~ lQ^^\t~'^'^ , where Kot is the rotating velocity in 
unit of cm s~^, t the age in seconds, and A a free parameter . 
Observations give A = 0.73 , according tolSkumanichI l^zj)) 
and A = 1.78, according to lSmithI l|l979f ). This estimates the 
rate of angular momentum loss by a low-mass rotating star 



in a binary system. Assuming that tidal coupling is per- 
fect, so that the wind-emitting star rotates synchronously 
with the orbital period, and neglecting the spin angular mo- 
mentum of each component relative to th e orbital angular 
mo mentum, we can t ake in to account the iTutukovl lll984) 
and llben fc Tutukovl (|l984^ formalism for the loss of angu- 
lar momentum by a MSW in the following form 



dlnJ 
At 
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(12) 



where J is the orbital angular momentum of the system 
in erg • s. Mi, a are the masses of both components in solar 
masses Mq, R2 is the radius of the donor in solar radii _Rq, 
and A is the semimajor axis of the system in solar radii Rq. 
So far we have only indirect evidence of the magnetic activity 
in W UMa-stars. Rapid period changes and disturbances in 
light curve s can be interpreted as being caused by big fiare s 
and spots jBinnendiikll977l: lMullanlll975l : IVan't VeeJl978ri . 



4 EVOLUTION INTO CONTACT 
4.1 The initial model 

Our initial model consists of two zero-age main-sequence 
(ZAMS) stars of Population I {X = 0.70, Z = 0.02) with a 
total mass of 1.8Mq (1.23 + O.57M0) and an initial mass 
ratio of 0.4634. The convective mixing length is twice the 
pressure scale hight. We take an initial period of 0.3766 
d and orbital angular momentum, J, of 5.161 x 10^^ erg-s 
[log(J/erg ■ s) = 51.713]. The initial separation between two 
components is about 1.8574 x lO" cm (i.e. 2.677?©). The 
initial model is a detached binary, and the surface of the 
primary lies only a short way inside its Roche lobe which it 
fills after about 7.3 x 10^ yr of nuclear evolution. 



4.2 Contact evolution without angular 
momentum loss 

4-2.1 Evolution without luminosity transfer due to mass 
transfer in semi-detached phase 

Our models in Paper I exhibit cyclic behaviour on a ther- 
mal timescale, with a period of about 10^ yr. Part of the 
time of each cycle is spent in a contact phase and part in 
a semi-detached phase. At first we only consider the lumi- 
nosity transfer by circulation in the common envelope under 
the assumption that no energy transfer due to mass transfer 
takes place in the semi-detached phase, the system evolves 
in thermal cycles without loss of contact. The changes of 
characteristic quantities (including the contact depth, sur- 
face temperature, central temperature, and central density 
of the primary and the secondary, together with the mass 
ratio and orbital period of the system) during the cycles are 
shown in Figure 1. 

As seen from Figure 1, beginning at the initial model, 
the primary fills its Roche lobe after about 7.3 x 10^ yr of 
nuclear evolution. It evolves into a semi-detached state and 
it loses mass to its companion. Since the addition mass to the 
secondary causes it to expand, and its effective temperature 
and luminosity to increase, after a total of about 1.5 x 10^ yr. 
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Figure 1. Without angular momentum loss and energy transfer due 
F, temperature, central temperature, central density of the primary 
mass ratio and orbital period of the binary] as a function of age. 

the secondary has swollen to fill its Roche lobe, so that the 
system evolves into a contact system. The contact binary 
formed at this point has masses of 1.0 and 0.8Mq, and a 
mass ratio of about 0.8. Thereafter, the model exhibits cyclic 
evolution on a thermal timescale (about 6 x 10® yr) without 
loss of contact. Part of time of each cycle is spent in good 
thermal contact, with EW light curves and part in poor 
thermal contact, with EB light curves. As seen from Figure 
lb, the ratio of the timescales characterizing the phases in 
which the oscillating systems exhibit W UMartype (EW) 
and (3 Lyrae-type (EB) light curves is tew/teb < 1- 



It is seen in Figures Ic and d that the central tempera- 
tures and central densities of both components of the binary 
exhibit cyclic evolution. This causes the cyclic evolution of 
the core luminosities of both components of the binary. 




J I I I I I I I I I I I I l: 




5x1 10* 

Age (yr) 

to mass transfer, the evolution of some quantities [contact degree 
(solid lines) and the secondary (dotted lines) , together with the 



4-2.2 Evolution with luminosity transfer due to mass 
transfer in semi-detached phase 

We now consider the evolution of the binary system from the 
initial model with luminosity transfer due to mass transfer 
included in the manner described in Sect. 2.2.1. The sys- 
tem also evolves in thermal cycles without loss of contact. 
The evolution of some characteristic quantities is shown in 
Figure 2. It is seen in Figure 2 that the system evolves 
into a contact binary from a detached one after a total of 
about 1.35 X 10^ yr of evolution from the main sequence 
and the contact binary formed at this time has masses of 
about 1.045 and O.755M0 with a mass ratio of about 0.72. 
So this model spent shorter time in the semi-detached phase 
than the model in which the luminosity transfer caused by 
mass transfer had been neglected in the semi-detached state. 
This is because the energy source at the secondary's atmo- 
sphere, provided by accreting mass from the primary in the 
semi-detached phase, can hasten the expansion of the sec- 
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Figure 2. With the energy transfer due to mass transfer, but without AML, the same quantities as Fig. 1 against Age. 



ondary and so shorten the time spent in the semi-detached 
evolution. 

It is seen in Figure 2b that the ratio the timescales 
(tew/teb) characterizing the phases in which the oscillating 
systems exhibits W UMa-type (EW) and /3 Lyrae-type (EB) 
is still less than 1.0. It is much smaller than its observational 
value (5 to 6) derived by iRahunen (1981) and even smalle r 
than a smallest value (13/5) derived bv iRucihskil i2002l) . 
This suggests that the real W UMa systems may suffer an- 
gular momentum loss. Therefore, we restrict the following 
discussions to the possible evolutionary consequences of the 
models for W UMa systems with angular momentum loss. 



4.3 Contact evolution virith angular momentum 
loss 

4-3.1 Evolution with angular momentum loss by GR 

Assuming that the angular momentum of W UMa-type sys- 
tems is lost only by gravitational wave radiation, we obtain 
a model for W UMa-type systems that exhibits cyclic evolu- 



tion without loss of contact as the models mentioned above. 
The contact degree and surface temperature of both compo- 
nents, together with the mass ratio and the orbital period of 
the system, are shown in Figures 3a-6a. It is seen in Figures 
3a-6a that the model exhibits cyclic evolution on a thermal 



timescale (with a period of about 6 x 
of contact. 



10 yr) without loss 



As seen from Figure 4a, the ratio of the timescales 
(t"ew/teb) is still smalle r than 1.0, th is is still much smaller 
than the value given bv I Vilhul il98l[) according to observa- 
tio ns, and even sm aller than a smallest value (13/5) found 
bv lRuciriskil (|20o3)- This suggests that, although the gravi- 
tational radiation has a significant influence on the evolution 
of W UMa-type systems and other close binaries, it still can 
not prevent the cyclic behaviour of the model for W UMa 
systems. iVilhu (,1381i) had concluded that if the timescale 
of angular momentum loss for W UMa stars is shorter than 
the secondary's thermal timescale (except at extreme mass 
ratios where the primary nuclear evolution dominates) the 
system will soon fill the outer critical surface and then prob- 
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Figure 3. Evolution of the contact degree of the binary with angular momentum loss by GR or by MSW of the secondary (with the 
different values of A). 



ably coalesces into a single star, if the timescale is longer 
than the primary nuclear timescale, the effect will be negli- 
gible. If the angular momentum loss rate falls between these 
limits, it makes an important contribution to the contact 
evolution. The rate of angular momentum loss associated 
with gravitational radiation is 10~^^ to 10"^° yr~^, corre- 
sponding to a timescale of about 10^° to 10^^ yr. It is clear 
that the gravitational radiation life-time is long compared 
with the nuclear timescale of the primary. Systemic angular 
momentum loss due gravitational radiation must accelerate 
the evolution toward more extreme mass ratios (see Figure 
5a) rather than lead to orbital collapse, so that the effect 
of gravitational radiation on the cyclic evolution is almost 
negligible. The most likely angular momentum loss mechar- 
nism for W UMa systems is a magnetic stellar wind. There 
is considerable observational and theoretical evidence that 
stars which possess convoctivo envelopes and, at the same 
time, rotation can amplify primordial magnetic fields, emit 



a stellar wind and lose angular momentum as a consequence 
of the interaction between the fields and the wind. 



4-3.2 Evolution with angular momentum loss by MSW 

Equation (12) gives a rate of angular momentum loss (via 
MSW) of a close binary system in which the mass is trans- 
ferred from the star 2 (with a mass of M2 and a radius of R2) 
to star 1. This suggests that the rate of angular momentum 
loss by MSW is directly related to the direction of the mass 
transfer. Assuming that the mass is always transferred from 
the secondary to the primary in W UMa systems, and tak- 
ing A=1.8, 1.2, or 0.8, which are in its observational range, 
we investigate the infiuence of the parameter A on the evo- 
lution of the models for W UMa systems. In these cases the 
evolution of contact degree and the surface temperature of 
both components, together with the mass ratio and the or- 
bital period of the system, are also plotted in Figure 3b-6b, 
c, and d, respectively. 
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Figure 4. Evolution of the surface temperature (Teff) of both components of the binary with angular momentum loss by GR or by 
MSW of the secondary (with the different values of A). 



It is seen in Figures 3-6 that our models still exhibit 
cyclic behaviour without loss of contact on the thermal 
timescale of the secondary. This is similar to the case of 
constant angular momentum with the exception that the 
cycles are slowly shifted towards more extreme mass ratios. 
And that the smaller the parameter A is taken to be (i.e. 
the higher the rate of angular momentum loss is assumed 
to be), the longer is the period of cyclic evolution, and the 
larger is the ratio of the timescales characterizing EW and 
EB phases (tew and teb). If we take A = 0.8, the ratio of 
the timescales (tew and teb) has risen u p to 2 to 3 whic h 
is consistent with the result obtained bv iRucinskil i2002f) . 
Therefore, the most likely angular momentum loss mecha- 
nism for W UMa systems is MSW, and the value of A is of 
about 0.8 for W UMa systems if the mass is transferred from 
the secondary to the primary in W UMa systems. Mean- 
while, as seen from Figures 3-6, the mean mass ratio in a 
cycle becomes smaller and smaller as the evolution of the 
system proceeds, in any case with angular momentum loss. 



the mean orbital period of the system becomes longer and 
longer. I.e. the system drifts to smaller and smaller mass ra- 
tios and to longer and longer periods along contact branches 
of the cyclic evolution. So W-type W UMa systems evolve 
into A-type W UMa systems with low total masses. Since 
the systems with P < 0.35 d are apparently not observed 
to have EB-type light curves, the fraction of time spent in 
the broken-contact phase or in poor thermal contact phase 
is uncertain and possibly small. This suggests that the poor 
thermal contact or semi-detached state probably appears 
very rarely (or never appears) during the evolution of the 
real W UMa systems. In effect, the loss of angular momen- 
tum due to MSW is used here to gently squeeze the systems 
and thus prolong the contact phase of the TRO model to 
such a degree that the broken-contact or poor thermal con- 
tact phase appears very rarely (or never appears). 

If the mass is transferred from the primary to the sec- 
ondary, the rate of angular momentum loss should be ex- 
pressed as 
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We take A = 1.8 which is the largest value obtained by the 
observations of single stars, and corresponds to the smallest 
rate of angular momentum loss due to MSW for single stars, 
and find that the model soon fills its outer critical Roche lobe 
after a total of about 2.5 x 10^ yr of nuclear evolution. The 
system loses its mass and angular momentum through the 
outer Lagrangian point. Then it probably coalesces into a 
single star. Since the direction of mass transfer may change 
with time, we take A = 1.8, 1.2, and 0.8, respectively, and 
calculate the rate of angular momentum loss via Eq. (12) 
when mass is transferred from the secondary to the primary, 
otherwise via the Eq. (13). In these cases, the models rapidly 
drift to extreme mass ratios, and they cannot be in good 
thermal contact when the mass ratio of the models becomes 
sufficiently small. Because, once the mass ratio of the system 
is very small, the mass is transferred from the primary to 
the secondary when the model is in poor thermal contact, 



then the rate of angular momentum due to MSW of the 
primary is very large, and the model evolves towards good 
thermal contact rapidly, but the direction of mass transfer 
has been reversed before the models are in good thermal 
contact, and the rate of angular momentum loss due to MSW 
of the secondary is too small to further force the model to 
evolve towards deep contact. Therefore, once the mass ratio 
of the models is sufficiently small the model is always in 
poor thermal contact even if we take the parameter, A = 0.7 
(see Figure 7) although the models exhibit thermal cyclic 
behavior with a short period and a small amplitude. 

In effect, W UMa systems are surrounded by a common 
envelope and the loss of angular momentum is probably not 
related to the direction of the mass transfer and may be 
simultaneously caused by MSW of both components. We 
take A = 3.8, and the loss rate of angular momentum of the 
model is taken as the sum of loss rates predicted by Eqs 
(12) and (13). The evolution of some quantities is shown 
in Figure 8. The thermal structure of the two stars of the 
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Figure 6. Evolution of the orbital period (P) of the binary with angular momentum loss by GR or by MSW of the secondary (with 
the different values of A). 

binary (A = 3.8) during a phase of contact evolution and This suggests that the magnetic braking in W UMa systems 
of semi-detached evolution is shown in Figure 9. It is seen is weaker than that in the similar single stars or components 
in Figure 8 that our model almost maintains an essentially of non-contact binaries. The observations have shown that 
constant depth of contact, rather than cycle in evolution, both in X-r avs l|C ruddace fc Dupre3ll984| IVilh^Jll984^ and 
and that the system steadily drifts to smaller and smaller in the radio iHughes^^McLeM]^^84 . ^^8a) the W UMa sys- 
mass ratios. This suggests that the model evolves steadily tems are definitely less active than are similar single stars 
towards a system with the extreme mass ratio (see Figure or components of non-contact binaries. The lower coronal 
8e). Meanwhile, if we assume a somewhat smaller value of activity in W UMa systems is probably caused by the in- 
A (i.e. a larger rate of angular momentum loss), the system ternal rotation of their components or by their convective 
soon reaches the outer critical Roche lobe. It then loses its zone which is thinner than indicated by their surface tem- 
mass and angular momentum through the outer lagrangian peratures (Ruciiisk i 1986). The relative weakness of coro- 
point, so that it probably coalesces into a rapidly rotating nae in W UMa systems directly tells us that the largest- 
single star. On the other hand, if the parameter A is as- scale magnetic structures are smaller in W UMa systems 
sumed to be somewhat larger, the system exhibits cyclic and contain less plasma in them than in other similar but 
evolution, and the subsequent evolution of the system is non-contact stars. Therefore, the magnetic-braking effects 
similar to the case of constant angular momentum with the in W UMa systems are weaker than in other similar but 
exception that the cycles are slowly shifted towards more non-contact stars. Meanwhile, we can see from Figure 8b 
extreme mass ratios. In this case, the parameter A is much that the surface temperature of the secondary is unstable, 
larger than the value of similar single or non-contact stars. so that its luminosity is unstable. This is consistent with the 
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properties of W-type W UMa systems with variable light 
curves. Moreover, we can see from Figure 8c that the radius 
of the primary is expanding rapidly while the radius of the 
secondary appears to be essentially constant although the 
model is drifting to smaller and smaller mass ratios. This is 
very similar to the result predicted by observations jWand 
Il995r . This suggests that W phenomenon is likely caused by 
the expansion of the primary (Paper I) because the expan- 
sion of the radius of the primary depletes a part of the core 
luminosity so that the temperature of the secondary exceeds 
that of the primary. 

As seen from Figure 9, the outermost layers of the com- 
mon envelopes of both components are so similar that the 
difference is slight. Meanwhile, we can see in Figure 9b that 
the profile of the specific entropy in the surface layers of the 
secondary is very different from the result given in Paper I. 
It clearly rises up to a higher value in the surface layers of 
the secondary because of the energy transfer due to mass 
transfer. Figure 12 shows the distribution of the specific en- 
tropy in both components as a function of the mass (Figure 



12a) and the radius (Figure 12b) during a phase of semi- 
detached evolution and of contact evolution with efficient 
energy transfer. The semi-detached model has the same to- 
tal mass and mass ratio as the contact one. The fiat parts in 
outer envelopes of the components indicate the convective 
regions. It is seen in Figure 12 that the convective enve- 
lope of the primary in the contact system is slightly thicker 
than that of the primary in the semi-detached one and the 
convective envelope of secondary in contact system is much 
thinner than that of the secondary in the semi-detached one 
(see Figure 12b), but the total mass contained in the con- 
vective envelopes of the conta ct system, which is directl y 
related to the dynamo action jHurlev. Tout fc PolslEooil . 
is much less than that contained in the envelopes of the 
semi-detached one (see Figure 12a). We conclude that this 
is caused by luminosity transfer between the two compo- 
nents. Because part of the luminosity in the envelope of the 
primary is transported to the envelope of the secondary in 
the contact state, the temperature of the envelope of the 
primary decreases and the temperature of the envelope of 
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the secondary increases. The decrease in the temperature of 
the envelope of the primary would then drive the radiative 
temperature gradient, Vr, to increase and thence lead part 
of the radiative zone to be convective so that the the convec- 
tive envelope of the primary becomes thicker in the contact 
state. The increase in the temperature of the envelope of 
the secondary would drive the radiative temperature gradi- 
ent, Vr, to decrease and thence lead part of the convection 
zone to be in radiative equilibrium. As a result, in contact 
systems, the mass in the convective envelope of the primary 
increases by about 0.01 Mq, and the mass in the convective 
envelope of the secondary decreases by about 0.08 Mq in 
comparison with in the semi-detached one (see Figure 12a). 
Since the decrease of the mass in the convective envelope of 
the secondary is much more than the increase of the mass 
in the convective envelope of the primary, the total convec- 
tive mass in the envelope of the contact system is less than 
that in the envelopes of the semi-detached one. As a result 
of the less total mass contained in the convective envelopes 



of the W UMa systems, because of energy transfer, W UMa 
systems show lower activity than non-contact stars. 

We also investigate the evolution of the models which 
lose angular momentum steadily at the different rates, and 
find that the system can keep in shallow contact (with a 
degree of about 20%, see Figure 10) by a hypothetical an- 
gular momentum loss « 1.6 10~^yr~^, on the average, 
corresponding to a timescale of abou t 6 10^ yr whic h is in 
agreement with the result obtained bv lRahune nl lll98j) . If we 
assume a somewhat larger rate, the system soon reaches the 
outer critical surface. This was confirmed with a steady loss 
of angular momentum with a rate that was slightly larger 
than the critical loss rate (see also Robertson & Eggleton, 
1977). After reaching the outer critical surface, the system 
loses its mass and its angular momentum through the outer 
Lagrangian point s, so that the system probabl y coalesces 
into a single star llWebbinklll976l : iRahunenlllQS J) . On other 
hand, if the assumed loss rate is smaller than the critical 
rate, the system exhibits cyclic evolution with alternating 
good contact and poor contact phases on the secondary's 
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thermal timescale. In this case the subsequent evolution is 
similar to the case of constant angular momentum with the 
exception that the cycles are slowly shifted towards more 
extreme mass ratios. 



4.4 The evolution in the period-colour diagram 

As shown bv lEggenI (|l96llll967i) . the observed W UMa sys- 
tems are located in a strip of the period-colour diagram, the 
two observed boundaries of the period-colour diagram are 
written by^ahlcr (20023) as 



1.51ogTeff - logP = 5.975.. .6.15, 



(14) 



where P is the period in days. This region is limited by the 
solid lines in Figure 11. Also shown is the massive compo- 
nent of our model [log(J/erg ■ s) — 51.713, and A = 3.8]. 
As seen from Figure 11, our model is in a good agreement 
with the period-colour relation for observed W UMa sys- 




Log P(days) 

Figure 11. The evolution of the massive component of the bi- 
nary in the period-colour diagram. 
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tems. Moreover, our model is at least similar to Rahunen's 
(1981) one. 



5 DISCUSSION AND SUMMARY 

In semi-detached phase, the mass is transferred from the 
primary to the secondary, and the energy (including the 
gravitational energy, thermal energy, and radiative energy) 
in the transferred mass is also transferred to the surface 
of the secondary. The effect of the accreting matter on the 
radii of the accretors is related to their envelopes being ra- 
diative or convective, and calculations involving mass ac- 
cretion on to a star with thick convect ive or fully convec- 
tive envelope JWhvte fc EggletonllTosil . have shown that 
the accreting star shrinks (under the assumption that the 
entropy of the accreting matter is the s ame as the surface 
entrop y of the original star). However, ISarna fc Fedoroval 
il989h thought that the newly accreting matter has sub- 
stantially higher entropy than the original star and so the 



photosphere of the secondary in close semi-detached binaries 
is heated by the infalling matter. Due to this process, the 
accreting secondary sh ould significantly increase in radius 
JPrialnik fc Livid Il985^ . That is to say, the energy source 
at secondary's atmosphere provided by the accreting matter 
from the primary hastens the expansion of the secondary 
and so shortens the time spent in the semi-detached evolu- 
tion. Therefore, the luminosity transfer due to mass transfer 
plays an important role in the origin and evolution of con- 
tact binaries. Meanwhile, if the energy transfer due to mass 
transfer is applied in the surface layer of the secondary as a 
surface condition of the secondary, we can see that this en- 
ergy source has a significant influence on the thermal struc- 
ture in the surface layers of the secondary (see figure 9b). As 
seen from figure 9b, the entropy profile in the surface layers 
of the secondary in the semi-detached phase rises. It is very 
different from the result in Paper I (see Paper I's Figure 6), 
and the difference arises from the luminosity transfer due to 
mass transfer. 

W UMa systems have been divided by iBinnendiikI 
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l)l97rfl into A-type and W-type systems according to 
whether the primary minimum in the light curve is a transit 
or an occultation. This suggests that the temperature of the 
more massive component is higher in A-type systems while 
the temperature of the less massive star is higher in W-type 
ones, i.e, the so-called W phenomenon. It is seen in Figure 8c 
that the radius of the primary is expanding in good contact 
while the secondary almost maintains an essentially constant 
radius in low-mass W UMa systems. Since the expansion of 
the primary depletes a part of the core luminosity, the tem- 
perature of the secondary exceeds that of the primary. This 
suggests that W-phenomenon is caused by expansion of the 
primary (see Figure 8b and c) . It is similar to the statistical 
result derived by .Wang, 11995.) from the observational data. 

Our model exhibits cyclic behavior without loss of con- 
tact, with a period of about 6 x lO'' yr if we do not consider 
angular momentum loss from the system. And the ratio of 
the timescales characterizing the phases in which the oscil- 
lating systems exhibit W UMa-type (EW) and (3 Lyrae-type 



(EB) light curves is tew/teb < 1-0. This is much smaller 
than its observational value (5 to 6, Rahunen 1981 and 2 
to 3, Rucinski 2002). Meanwhile, all constant angular mo- 
mentum loss models need a formation mechanism which can 
produce unequal components because, with constant angu- 
lar moment um, it is not po ssible to co ver a large range of 
mass ratios llVilhulll 981: Rah unenlllQSlTl . This suggests that 
the real W UMa systems should suffer angular momentum 
loss during their evolution. The rate of angular momentum 
loss associated with gravitational ra diation is about 10~^^ 
to 10~^° yr~^ for W UM a systems (iRobertson fc EggletonI 
I1977I: IWebb'in3ll975l iRa h unen 1^. This is two orders of 
magnitude smaller than the rate of the hypothetical angular 
momentum loss needed to keep the systems in good ther- 
mal contact. The effect of gravitational wave radiation on 
the cyclic evolution is so small that the model with angu- 
lar momentum loss due to GR is not consistent with the 
observations of W UMa systems. The most likely angular 
momentum loss mechanism is MSW. In fact, there are some 
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indications of_strong magnetic activity in W UMa systems. 
IVan't Veeil il979h thinks that magnetic braking may con- 
trol the evolution of solar-type contact binaries. Eqs. (12) 
and (13) give the rate of angular momentum loss by MSW 
of the secondary and the primary, respectively, and the ob- 
servations give the value of parameter, A = 0.73 to 1.78 for 
single stars. According to Eqs. (12) and (13) the rate of an- 
gular momentum loss caused by MSW is directly related to 
the direction of the mass transfer in a close binary system. 
If we assume that the mass is always transferred from the 
secondary to the primary, and take A = 0.8, 1.2, and 1.8, re- 
spectively. We find that our models exhibit cyclic evolution 
which is similar to the case of constant angular momen- 
tum with the exception that the cycles are slowly shifted 
to the more extreme mass ratio, and the larger the rate of 
angular momentum loss is assumed to be, the more rapidly 
the model evolves towards the extreme mass ratios, and the 
larger is the ratio of the timescales (tbw/teb). If the param- 
eter, A, is taken as 0.8, the ratio of the timescales (tew/teb) 
is abo ut 2 to 3 which is in agreement with a observational 
value llRucinskill26o3) . Meanwhile, we can see from Figures 
3-6 that the mean mass ratio of a cycle becomes smaller and 
smaller and the mean period becomes longer and longer as 
the evolution of the system proceeds. This indicates that it 
slowly evolves towards a system with an extreme mass ratio 
and a long period. If the mass is assumed to be transferred 
from the primary to the secondary, even if the parameter, A, 
is taken to be 1.8, the model soon fills its outer Roche lobe. 
The star would then lose its mass and angular momentum 
through the outer Lagrangian point, and the system would 
soon evolve into a single star. If the direction of the mass 
transfer is changed with time, and the angular momentum 
loss of the system is caused by MSW of the secondary when 
the mass is transferred from the secondary to the primary or 
by MSW of the primary when the mass is transferred from 
the primary to the secondary, no matter how the parameter 
A is chosen between 0.73 and 1.78, the model cannot be in 
good contact when it evolves to a system with a sufficiently 
small mass ratio. 

Our models for low-mass W UMa systems show that 
the magnetic braking effects in W UMa systems are likely to 
be weaker than that in other similar but non-contact stars. 
This is in agreement with observations. Althou gh there 
were both observationa l and theoretical indications llRobertj 
Il974l: iKnobloch eraP Il982i) of a power law {B = BqQ.^ , 
where fi is the Keplerian angular velocity, B the average 
poloidal field, with the most probable value of 5 ~ 1). 
However, observation has shown that the strengths of the 
chromospheric and transition-region emissions increase very 
slowly with the increase of the inverse to Rossey number t/P 
jRucmski 19861: IStepieri 1995ft . This can be explained as a re- 
sult of the total saturation of stellar surfaces by small active 
regions. Meanwhile, the observations have shown that, both 
in X-rays and in the radio, the W UMa systems are definitely 
less active than a re similar single stars or components of non- 
contact binaries llRuciiiskilll986ll . The lower coronal activity 
in W UMa systems is likely caused by the internal rotation of 
their components or because their convective zone is thinner 
than indicated by their surface temperatures. The observa- 
tions directly tell us that in W UMa systems the largest-scale 
magnetic structures are smaller and contain less plasma in 
comparison with the similar but non-contact stars so that 



the efficiency of the angular momentum loss in W UMa sys- 
tems may be lowered in comparison with non-contact stars 
of the same t/P. Figure 12 shows the distribution of the 
specific entropy of both components during a phase of a 
semi-detached state and of a contact state. It is seen in Fig- 
ure 12 that the convective envelope of the primary in the 
contact system is slightly thicker than that in the similar 
semi-detached one and the convective envelope of the sec- 
ondary in the contact state is much thinner than that in 
the semi-detached state, but the total mass contained in the 
convective enve lopes of contact systems, which is related to 
dynamo action (iHurlev. Tout & Pols 2002), is less than that 
contained in the convective envelopes of the semi-detached 
stars because of energy transport from the primary to the 
secondary (see Figure 12). This is in a good agreement with 
the prediction of observations mentioned above. As result of 
the less total mass in convective envelopes of W UMa sys- 
tems, W UMa systems show the lower activity in comparison 
with the single or non-contact stars. 

Our models indicate that W UMa systems evolve into 
contact binaries with the extreme mass ratios or even into 
single stars i f they suffer angular momentum loss due to 
MSW. IVilhul Jl98lft thought that W UMa systems most 
probably evolve into single stars, but the way (and the 
timescale) in which this evolution proceeds is not clear. How- 
ever, it is argued that due to a longer timescale to reach a 
contact, and also to the longer timescale of contact evo- 
lution, the ultimate appearance of low-mass binaries after 
merging of both components will be distinctly different from 
that of massive binaries. After an initial mass exchange of 
a massive contact system, it is very likely that the main- 
sequence (MS) lifetime of a more massive component be- 
comes shorter than the duration of the contact phase of evo- 
lution, and that the coalescence results in the formation of 
a fast-rotating giant, similar to FK Com. In the case of low- 
mass binary, the MS lifetime of a more massive component 
may still be longer than the duration of the contact phase. 
That would result in the formation of a single fast-ro tating 
MS star, similar to single blue straggler llStepienlll995ft . The 
contact systems with most extreme mass ratios (few known 
due to the low probability of eclipses) are important for test- 
ing models of contact binary evolution. But most W UMa 
stars with low mass ratios (such as RR Cen and e CrA, 
see Vilhu 1981) seem to have an evolved primary and both 
components should deviate from thermal equilibrium. 

Although we can keep our models in good contact un- 
der some assumptions, there is still one observational point 
which seems to be in contradiction with the scenarios pre- 
sented above, because it produces a peculiar mass ratio dis- 
tribution showing a large excess of the systems with small 
mass ratios. This is due to the fact that the mass trans- 
fer occurs on the thermal timescale of the secondary, which 
grows strongly as the secondary mass and therefore as the 
mass ratio decreases. Thus systems with large meiss ra- 
tios initially evolve in less than 10* yr towards mass ra- 
tios less than about 0.4 (see Figures 8 and 10). The sce- 
narios in which contact-binary evolution continues into a 
single star on the secondary thermal timescale produce a 
rather steep period-distribution. It can be shown quite eas- 
ily that the resulting mass ratio distribution does not agree 
with the mass ratio distribution deduced from observa- 
tion iVan't Veeilll978l : IVilhulll98ll) . Meanwhile, W UMa- 
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type contact binaries are found in stellar aggregates of a 
widely differing age, in very old clusters like NGC188 and 
in moderately old ones like M67 and Praesepe. This cannot 
be interpreted by these models. The decrease in the coro- 
nal (X-rays and ra dio) activity of contact binaries, which 
is indeed observed llCniddace fc DuDreelll984t IVilhxJll984l: 
iHughes fc McLeanlllOslT is used as an argument that the 
AML efficiency in contact is relatively low and the contact 
stage is considerably prolonged relative to adjacent stages. 
This small modification to the AML models is capable of 
explaining why so many different contact binaries are ob- 
served in old systems like NGCf88. Therefore, a fee dback 
mechanism is crucial to the success of these models. IVilhrJ 
il98lfl envisaged a situation in which increasing depth of 
contact causes increased mixing in the common envelope 
which tends to bury the strong surface magnetic field. The 
magnetic braking is then weakened, leading to a decreased 
loss rate of angular momentum and decreased depth of con- 
tact again. It is clear that this picture is very crude and 
speculative. On the observational side one should try to mea- 
sure the magnetic fields and to study the chromospheres and 
coronae more closely in the UV- and X-ray region. On the 
theoretical side one should clarify the differential rotation 
and dynamo action in rapidly rotating stars in connection 
with the circulatio n pattern in t he common envel ope, a task 
which is not easy llVilhulll98ll) . lRuciiiskil Jl982l) explained 
the period gap at periods just longer than those for con- 
tact systems by a phase of very rapid magnetic braking in 
detached binaries which ends when the stars come into con- 
tact. However, it is unclear that this mechanism is sensitive 
enough to maintain an essentially constant depth of contact. 



The spin angular momentum of both components and 
the loss of the system's mass have not been considered in 
our present work. Indeed the spin angular momentum of 
both components is not small compared with the orbital 
angular momentu m when the mass ratio becomes extreme 
jVan't Veeilll979D . If the spin angular momentum of both 
components is considered, angular momentum loss required 
to spin up the stars to keep them in corotation as the orbit 
shrinks, so the spin angular momentum of both components 
increases and the orbital angular momentum of the system 
decreases, the system might then suffer Darwin's instabil- 
ity (i.e., the tendency for an orbit to desynchronize if the 
spin angular momentum of both stars is more than a third 
of the orbital angular momentum, Eggleton & Kiseleva- 
Eggleton 2001) when so much angular momentum has been 
removed that there is no longer sufficient in the orbit to bring 
the stars into corotation and evolution proceeds on a tidal 
timescale (i.e., rtid = , where the dot indicates a time 

"spin 

derivative, Ogpin is the spin angular velocity, Bagot 1996). 
As a result, angular momentum loss would ultimately lead 
to a lack of corotation and thence to a failure of the Roche 
model. Meanwhile, magnetic braking is always accompanied 
by the mass loss from the system, and the mass loss has a 
significant influence on the evolution of W UMa systems. 
We will consider the spin angular momentum and mass loss 
of the system to construct a true non-conservative model for 
low-mass W UMa systems in our future work. 
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